Neel to Spin-Glass-like Phase Transition versus Dilution in Geometrically Frustrated 

ZnCr2_2xGa2x04. 
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ZnCr204 undergoes a first order spin-Peierls-like phase transition at 12.5 K from a cubic spin 
liquid phase to a tetragonal Neel state, [ll] Using powder diffraction and single crystal polarized 
neutron scattering, we determined the complex spin structure of the Neel phase. This phase 
consisted of several magnetic domains with different characteristic wave vectors. This indicates 
that the tetragonal phase of ZnCr2_2xGa2x04 is very close to a critical point surrounded by many 
different Neel states. We have also studied, using elastic and inelastic neutron scattering techniques, 
the effect of nonmagnetic dilution on magnetic correlations in ZnCr2_2xGa2x04 (x=0.05 and 0.3). 
For x=0.05, the magnetic correlations do not change qualitatively from those in the pure material, 
except that the phase transition becomes second order. For x= 0.3, the spin-spin correlations become 
short range. Interestingly, the spatial correlations of the frozen spins in the x=0.3 material are the 
same as those of the fluctuating moments in the pure and the weakly diluted materials. 
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I. INTRODUCTION 

There has been a long standing fascination in the 
physics community with placing antiferromagnetically 
coupled spins on lattices with triangular motifs. 0, Q 
In two dimensions, one can consider the simple trian- 
gular lattice. When classical XY spins are placed on 
this lattice, it orders at zero temperature. For many 
years, the excitations above this ground state provided 
a playground for studying topological excitations. These 
chiral excitations destroy the ground state at finite tem- 
peratures and experimental investigations of these sys- 
tems are still underway. One can increase the com- 
plexity of this problem by considering the case of the 
even less constrained system of corner sharing trian- 
gles, the kagome latticed, H, i, [13, HH There have 
been many theoretical questions about the nature of the 
ground state [l^ [l^, [3, EB] even for classical spins on 
this lattice. [^,^^1 13| Experimentally, good experimen- 
tal realizations of this system have been hard to ob- 
tain. Recently, however, single crystals of Fe-jarosite 
that realize the kagome lattice with classical (S = 5/2) 
spins were grown and detailed neutron scattering studies 
have been performed. [l6| More recently, ZnCu3(OH)6Cl2 
was found as a good model system for the quantum 
kagome antiferromagnet. 17, 18, JJ,, 20., 211 Unfortuately, 
it turned out that there is 5-10% site switching of Zn^+ 
and Cu^^ ions.[2l| which obscures the quantum physics 
of the kagome antiferromagnet. 

In three dimensions, when spins form a network of 
corner sharing tetrahedra, it leads to a macroscopi- 
cally degenerat e g round state for classical as well as 
quantum spins. [12 [23| Theoretically novel low temper- 
ature properties are expected to appear. For exam- 
ple, quantum spin liquid phases, fractional excitations, 
or non-Ginzburg-Landau phase transitions. Experimen- 
tally, spinels AB2O4 have attracted lots of attention be- 



cause the B sublattice forms a network of corner-sharing 
tetrahedra. In the spinel, the B site cations are octahe- 
drally coordinated by six oxygens and neighboring BOg 
octahedra share an edge. Thus, when the B site is occu- 
pied by a transition metal ion with t2g electrons, the sys- 
tem can realize the simple and most frustrating Heisen- 
burg spin Hamiltonian, H = J ^ • Sj with uniform 
nearest neighbor interactions. 

ACr204 (A = Zn[H, Cd^, HgH^) realizes the most 
frustrating Hamiltonian because the t2g orbital of the 
Cr'^^ {3d^) is half filled and the nearest neighbor inter- 
actions due to the direct overlap of the neighboring t2g 
orbitals are dominant and spatially uniform. 27] In com- 
parison, in the case of AV2O4 where the V^+(3rf^) ion has 
an orbital degeneracy, a Jahn- Teller distortion can occur 
at low temperatures, which makes the vanadates effec- 
tively one-dimensional spin chain systems. [28l. [29l. [sol [Blj 
Several novel discoveries have been made in ACr204. 
For instance, collective excitations of local antiferromag- 
netic hexagonal spins were found in the spin liquid phase 
of ZnCr204 that embody the zero-energy excitations 
amongst the degenerate ground states. [24, 32] Unfortu- 
nately, the lattice of ACr204 is not infinitely firm and 
it distorts at low temperatures to lift the magnetic frus- 
tration. The novel three-dimensional spin-Peierls phase 
transition, i.e., the lattice instability driven by mag- 
netic interactions, occurs and drives the system into a 
Neel state. The lattice distortion can occur in different 
forms, depending on details of the crystal environment: 
tetragonal /4m2 symmetry for ZnCr904[33j. tetragonal 
lAi/amd for CdCr204[33|, and orthorhombic Fddd for 
HgCr204[2B, [S^]- When an external magnetic field is ap- 
plied to the Neel state, the half-magnetization plateau 
states appear in CdCr204[3^ and HgCr^04[26l. due 
to the field-induced lattice instability 3^ [13] ■ 

In this paper, we investigated the nature of the 
3D spin-Peierls transition in the chromite by perform- 
ing elastic and inelastic neutron scattering measure- 
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ments. Neel phase survives up to x ^ 0.2. For 0.2 < 
X < Xc with 1 — Xc — 0.390(3) being the percolation 



FIG. 1: Octahedral B-sites of spinel AB2O4 form a network 
of corner-sharing tetrahedra. 



ments on ZnCr2-2xGa'2x04 for x=0, 0.05, and 0.3. 
Our principal results are the following. For x=0, the 
Neel state has four characteristic wave vectors, k — 
(l,0,0),(i,i,i),(l,0,i) and (i,i,0).[33 The large size 
of the magnetic unit cell (64 Cr'^+ ions) has made it dif- 
ficult to uniquely determine the spin structure of this 
system. We have determined the spin structure, em- 
ploying powder diffraction, single crystal polarized neu- 
tron diffraction data and a systematic group theoretical 
approach [53. We find that the system is composed of 
three types of domains whose relative fractions vary from 
sample to sample. The dominant domain is a multi- 
k structure with k — (1,0, i) and (^,^,0). This spin 
structure is coplanar and noncoUinear with spins pointing 
along either the a or 6 axis with each tetraheron having 
two pairs of antiparallel spins to have zero net moment. 
The k = domain has a rather simple spin struc- 

ture. The spinel lattice can be decomposed into alternat- 
ing kagome and triangular layers when viewed along the 
< 111 > direction. In this spin structure, the spins in 
the kagome layer order in the "q=0" configuration. The 
spins in the triangular layer point along the < 111 > 
direction and are parallel within a layer. Spins in alter- 
nating layers are antiparallel. The k = (1,0,0) domain 
has a collinear spin structure with spins parallel to the z- 
axis, as in ZnV204.[29] In each tetrahedron the net spin 
is zero. 

The effect of site disorder on the magnetic corre- 
lations and phase transition in ZnCr204 by doping 
nonmagnetic Ga ions into Cr sites has been previ- 
ously studied by Fiorani et al. using bulk property 
measurements and neutron powder diffractions in the 
ZnCr2-2xGa2x04 series [40|. Fig. 2 shows the phase 
diagram that they have constructed from the measure- 



threshold for the corner-sharing tetrahedra [4l|, the sys- 
tem exhibits spin-giass-like properties in bulk susceptibil- 
ity measurements. However, the nonlinear susceptibility 
of ZnCr2-2xGa2x04 (x=0.2) does not display the diver- 
gence expected of an ordinary spin glass [43|. This sug- 
gests that the low temperature phase is not an ordinary 
spin glass. 
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FIG. 2: Phase diagram for ZnCr2-2xGa2x04 . Neel tem- 
perature Tjv and spin freezing temperature T/ represented 
by circles are the data obtained by bulk susceptibility mea- 
surements reported in Ref. [4^. Tn represented by squares 
and in the inset are obtained by our neutron scattering 
measurements, which is discussed in the Sections IV and V. 



Our data show that in the weakly diluted 
ZnCr2-2xGa2x04 (x=0.05) Neel ordering occurs 
with the same spin structure as the parent compound. 
However, the ordering now develops gradually and the 
phase transition becomes second order. The appearance 
of the Neel ordering and the cubic-to-tetragonal struc- 
tural phase now also proceed in a second order manner. 
This consanguinity of the order of the structural phase 
transition and the appearance of long-range magnetic 
ordering supports our interpretation that the phase 
transition is magnetically driven. 

For ZnCr2_2xGa2x04 (x=0.3), the magnetic long range 
order is replaced by static short range order even though 
the Cr concentration, 1 — x = 0.7, is well above the per- 
colation threshold, 0.390. Interestingly, the spatial cor- 
relations of the frozen spins in the spin-glass sample are 
the same as those of the fluctuating moments present in 
the pure and weakly diluted materials. Magnetic neu- 
tron scattering intensity, I{Q), vanishes as Q — *■ and 
has a broad peak at Qc — 1.5 with fuU-width-of- 
half-maximum (FWHM) of k = 0.48(5) A'^ This in- 
dicates that fundamental spin degree of freedom in the 
corner-sharing tetrahedra involves an antiferromagnetic 
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hexagonal spin loop with zero net moment, which distin- 
guish the geometrically frustrated magnet from an ordi- 
nary spin glass. 

The structure of this paper is the following: In section 
II, we describe the experimental details of material syn- 
thesis and the neutron scattering techniques that were 
used. In Section III, we explain the determination of 
the spin structure of ZnCr204. In Section IV, we dis- 
cuss inelastic neutron scattering data on the material and 
how the spin freezing and short range spin correlations 
in the diluted compound resemble those in other frus- 
trated magnets. This paper concludes with a discussion 
and summary in section V. 



II. EXPERIMENTAL DETAILS 

Three 20 g powder samples of ZnCr204 , one 
ZnCri gGao.i04 sample and ZnCri 4Gao.604 sam- 
ple were prepared by the standard solid state reaction 
method with stoichiometric amounts of Cr203, Ga203 
and ZnO in air. Neutron powder diffraction measure- 
ments performed on the samples at the National Insti- 
tute of Standards and Technology (NIST) BTl diffrac- 
tometer show that the samples were stoichiometric single 
phase spinels with the exception of one ZnCr204 sample 
(sample 2) which had a minority phase of 1% f.u. un- 
reacted Cr203. The results of the structural refinement 
are summarized in Table I. The ZnCr204 samples will 
be denoted by sample 1, sample 2, and sample 3 in this 
paper. 

A 0.1 g single crystal of ZnCr204 was grown by the 
chemical transport method and used for polarized neu- 
tron diffraction measurements at the NIST cold neutron 
triple-axis spectrometer, SPINS. Spectrometer configu- 
ration was guide-PG(002)-Be-Pol.-40'-Samp.-Flip.-Pol.- 
40 -PG(002)-Det. The sample was mounted such that 
the scattering plane were the (hkO) and (hOl) zones due 
to twinning. A vertical guide field was applied. The 
polarization efhciency was determined by measuring the 
scattering intensities of a nuclear (2,2,0) Bragg peak with 
the flipper on and off. Correction for the finite polarizing 
efhciency, 0.85(1), was made|4J]. 

For inelastic neutron scattering measurements on pow- 
der samples, we utilized a multiplexing detection system 
of SPINS consisting of a flat analyzer and a position- 
sensitive detector. The details of the experimental setup 
are reported elsewhere [l|. High angle backgrounds were 
measured by defocusing the analyzer while low angle 
backgrounds from air scattering were measured by ex- 
tracting the sample from the cryostat. The absolute ef- 
ficiency of the detection system was measured using in- 
coherent elastic scattering from vanadium and nuclear 
Bragg peaks from the samples. The corresponding cor- 
rection factor was applied to the background subtracted 
data to obtain normalized measurements of the magnetic 
scattering cross section [431]. 



III. ANTIFERROMAGNETIC LONG-RANGE 
ORDER IN ZNCR2O4 



A. Powder diffraction data 
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FIG. 3: Powder diffraction data from ZuCt^Oa at 25 K 
and 2 K taken on the powder diffractometer BT 1 at NIST. 
The line through the data in (a) shows the Rietveld fit to 
the crystal structure of which parameters are shown in Ta- 
ble 1. The line in (b) is the fit to the crystal structure 
and the magnetic structure that are explained in the text 
for sample 1. The upward pointing arrows with plain line 
are magnetic reflections which belong to k = (i, i, 0). From 

the left, they are (i, i, 0), (i, i, 1), (|, |, 0), (|, |, 1), (i, i, 2) 

and (|, 1,0), (i, 1,2) and (i,f,0), and (i,f,l). The 
upward arrows with dashed line are magnetic reflections 
which belong to k = (1,0, |). From the left, they are 
(1,0, i), (1,0, I), (1,2, i), and (2, 1, |) and (1,0, f). The 
downward pointing plain arrows are from the k = (i, |, i) 
family of magnetic reflections. From the left, they are 

(iif). (I'l'D' ™d (|,f,^). The downward pointing 
dashed arrows are from the k = (1,0,0) family of magnetic 
reflections. From the left, they are (1,1,0), (2,1,0), and 
(2, 1, 1). Figures (c) and (d) show the data and fits for samples 
(2) and (3) respectively, also discussed in the text. 
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Fig. 3 (a) shows T = 25K > Tjv diffraction data from 
ZnCr204 with the Rietveld fit superimposed. In addition 
to the nuclear Bragg reflections, there is a broad peak 
centered at Q ~ 1.5 A~^. This broad peak is due to dy- 
namic spin fluctuations and will be discussed in the Sec- 
tion IV. Below Tc, the broad peak weakens and magnetic 
Bragg reflections appear, indicating a long-range mag- 
netic ordering. These diffraction patterns are consistent 
with those observed previously from ZnFen iCri gO^fist. 
Indexing these magnetic reflections indicates that the 
magnetic unit cell consists of four chemical formula units 
(64 magnetic Cr^+ ions) which can be characterized 
by four wave vectors, k = (i,i,0), (1,0, i), (^,^,^), 
(1,0,0) It is impossible to uniquely determine spin struc- 
ture for a system with such a large magnetic unit cell only 
from its powder diffraction pattern. Oles proposed an in- 
plane spin structure for ZnFen iCri QOzt|45l| and Shaked 
et al. a non-inplane structure for MgCr204[46'|. Appar- 
ently, as we will show in Section III. C, there are numer- 
ous spin structures that can explain the neutron powder 
diffraction data equally well. To obtain more restrictive 
information for the spin structure, we have performed po- 
larized neutron diffraction on a single crystal of ZnCr204. 
Our polarization study focuses on the k = (i, ^,0) and 
(1, 0, ^) family of magnetic reflections. 



B. Polarized neutron diffraction data from a single 
crystal 



from the incomplete instrumental polarization efRciency 
of 0.85(1). For magnetic Bragg reflections that belong to 
the k = (i, i,0) family, Eq. (1) becomes 
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FIG. 4: Rocking scan through (a) a nuclear Bragg reflection 
at (2,2,0) and (b) a magnetic Bragg reflection at (i,— 1,0) 
obtained from a single crystal of ZnCr204 at T = 1.7 K. Open 
circles are the NSF data and filled ones are the SF data. 



This material undergoes a cubic-to-tetragonal struc- 
tural phase transition with c < a at T^v = 12.5 K[i[. 
Because of tetragonal twinning, below T/v a wave vec- 
tor transfer Q = {Qx,Qy,Qz) in the laboratory coordi- 
nate system represents {Qx/a* ,Qy/a* ,Qz/c*) = {h,k,l) 
and {Qx/a*,Qy/c*,Qz/a*) = {h,l,k) in different crystal 
twin domains. In the configuration with a vertical guide 
field, the non-spin- flip (NSF) and spin-flip (SF) scatter- 
ing cross sections, aNSF and cfsf become pTl |48| 



(^NSF — (^N 
C^SF 



a 



MA 



'MA 
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Here ctat is structural scattering cross section, and aM is 
the magnetic scattering cross section, (Jm oc (1 — Q ■ 

S)\Fm{Q)\'^ ■ In <^NSF, we neglected the interference 
term between the nuclear and the magnetic scattering 
amplitude because the reflections considered here are ei- 
ther purely nuclear or purely magnetic. Fm(Q) is the 
magnetic structure factor Fm{Q) = ^gF{Q)T,d < Sd > 
gjQ d T^\yQj-Q F{t) is the magnetic form factor of Cr'^+[5l[. 
The perpendicular sign in aM± is to emphasize that only 
the spin components perpendicular to the wave vector 
transfer Q contribute to scattering. Fig. 3 (a) shows 
the NSF and SF scattering intensities obtained at a nu- 
clear (2,2,0) Bragg reflection. Only NSF scattering is 
expected for this nuclear Bragg reflection and the con- 
tribution in the SF channel due to the contamination 



As shown in Fig. 4 (b), The magnetic (^,|,0) 
Bragg reflection has dominantly SF intensity and a 
weak signal in the NSF channel. The weak NSF inten- 
sity is due to contamination from incomplete polariza- 
tion. We investigated five magnetic Bragg reflections: 
(i,i,0),(i_,|,0),(f,f,0),(i,|,0), and (|,f,0). After 
the correction for the incomplete polarization, at these 
reflections only SF scattering is present. This means that 
the spins are in the ab-plane, a%j — 0. 

For the magnetic Bragg reflections that belong to 
k — (1,0,^), there are NSF as well as SF scattering. 
The ratio of SF to NSF scattering intensity, ctsf I o'nsf^ 
is different at different Q, as shown in Fig. 5. Table I 
lists the ratios for three different reflections. The ratio, 
c^sF/crNSF, increases as I increases. This information 
places a strict restriction on any possible spin configura- 
tion for ZnCr204. For these reflections, Eq. (1) becomes 



<7nsf — c 
asF = c 
Therefore, using a%j = 0, 
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The fact that the measured (Jsf/(^nsf follows 1 — 
within the experimental accuracy indicates = 
for the k = (1,0, i) reflections. 
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functions V'x^ for k 
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FIG. 5: Rocking scan through magnetic Bragg reflections (a) 
at (1, ^,0) and (b) at (1, §,0) obtained from a single crystal 
of ZnCr204 at T = 1.7 K. Open circles are the NSF data 
and filled ones are the SF data. 



C. Group theoretical approach to determination of 
spin structure 

We have employed a group theoretical approach devel- 
oped by Izyumov et al. |50i] to determine the spin struc- 
ture. The basic idea of the method is that any magnetic 
structure with a characteristic wave vector k can be ex- 
panded in terms of basis functions, ip^^ j of irreducible 
representations of the spin space of the crystal Gk which 
is a subgroup of the crystal space group g[50I|. 



c{k} 



(5) 



Here Oj represents a magnetic ion at site j in zeroth prim- 
itive cell. L runs over the arms of the star k^ and A 
over irreducible representations of the star arm k^ . The 
star of a wave vector k, {k}, is the set of nonequiva- 
lent vectors that can be obtained by acting on k with 
an element of the crystal space group 5 G G. For in- 
stance, the star {k} = {(1,0, i)} has six arms: k^ — 

(l,0,i),(l,0,i),(l,i,0),(l,i,0),(i,0,l), and (i,0,l). 
Once Sq^^ is determined, all spins at other primitive cells, 
Snj, can be derived bvfsoj 



■^"i ^ E exp(ik • tn)S^ 



(6) 



where t„ is the translation vector for the n cell from the 
zeroth primitive cell. How to obtain the basis functions 
of the irreducible representations, V'a^ for a given Gk 
has been explained in a great detail in a book by Izyu- 
mov et aZ.[50| and will not be repeated here. The basis 



listed in Table II and III respectively. Note that ip^^ 
are complex but their simple superposition at the two 
arms Ciip^^'^ + C2tp^^'^ with k2 = — ki can generate a 
real function, provided that the coefficients Gi and G2 
are appropriately selected. Table IV and V list such su- 
perpositions for k = (i,i,0) and k = (1,0,^), respec- 
tively. The superpositions of two irreducible represta- 
tions ,Gi'0''^^ + C2il'^^'^ , do not yield nonzero spins for 
all 16 magnetic ions in a chemical unit cell. Instead the 
superposition yields four nonzero spins for k — (i, ^,0) 
and eight nonzero spins for k = (1,0, i) in a chemical 
unit cell. This means that to put all 16 nonzero spins 
into a chemical unit cell we have to consider at least four 
of the Ciij^'^^+C2iJ^^^ listed in Table III for k = (i, i, 0) 
and at least two of Ciip^^^ + C2ip^^'' listed in Table IV 
for k = (1,0, ^). Obviously the number of such combi- 
nations is very large. To narrow down the possible spin 
structure, we used the constraints that were obtained 
from polarized neutron diffraction data: (1) spins are 
in the ab-plane, Sc = 0, and (2) = \F^^\ for the 

(1,0, ^) family reflections. We also assumed that (3) all 
spins have the same magnitude and that (4) all tetrahe- 
dra have zero net spin. 

Since for this domain, there are two characteristic wave 
vectors involved in the Neel state of ZnCr204, we rewrite 



Eq. (5) to separate S", 



{k} 



S, 



{k} 



into two components; 

ki(l,0,i) , „ki,(i,i,0) 



Oj 



„{k(l,0,i)} „ 



L 



(7) 



■,o)} 



Note that S'{k(i,o,i)} contribute only to the k = (1,0, i) 

family reflections and S^^^^'^''^''^ only to the k = {^,^,0) 
reflections. Therefore we can obtain the two components 
separately. 

First, let us consider the k = (1,0, i) reflections. We 
examined all possible combinations of the superpositions 
listed in Table V and found 24 different spin configura- 
tions which can be divided into three categories shown 
in Fig. 5. Fig. 6 (a) shows a non-coUinear spin config- 
uration in which spins are along either a-axis or b-axis. 
Fig. 6 (b) shows a non-coUinear spin configuration with 
spins along (1,1,0) or (1,1,0), and Fig. 6 (c) shows a 
coUinear spin configuration along (1,1,0). A collinear 
spin configuration along a-axis or b-axis is ruled out by 
the constraint = \F^j\. 

Fig. 8 shows four prototypes of k = (i, i, 0) spin con- 
figurations in which all tetrahedra satisfy the antiferro- 
magnetic constraint to have zero net moment: (100) type 
collinear and noncoUinear spin configurations and (110) 
type collinear and noncoUinear spin configurations. 

Now, it is possible that ZnCr204 has two magnetic 
domains: one with k = (1,0,^) and the other with 
k = (i,i,0). However, to explain the neutron pow- 
der diffraction data, the population of the two domains 
as well as the ordered moment have to be exactly the 
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FIG. 6: Prototypes of possible spin structures in a chemi- 
cal unit cell for k = (1,0, |) which satisfy the conditions de- 
scribed in the text: (a) V''^^" -FV''^^'^ +itp^^'''ii -itp^^''i2 , (b) 

{l + i)i>^^^^^ + (1 - i)t/'''2^i" + (-1 + + (-1 - , 

and (c) +i>^^^^^ +iTP^^^'ii -ii>^^^'i2 -iip^^^^i -ji/>''2^22 _ 

^kir^i _^^k24^ Here ki = (1,0, i) and k2 = (1,0, i). Mag- 
netic unit cell is doubled along the c-axis and spins change the 
sign in the chemical unit cell displaced by (0,0,1). A shaded 
square represents a tetrahedron formed by four Cr'^+ ions. 
Symbols represent z-coordinates of the magnetic Cr^"*" ions. 



same for both domains. Furthermore, another spinel 
ZnFe204 also magnetically orders at low temperatures 
and the Neel state has only a single characteristic wave 
vector k = (1,0, i). We believe it is more likely that 
the two characteristic wavevectors participate in the or- 
dering of all the spins in ZnCr204. Then, the resulting 
spin structure would be a summation of 5'^'^'-^'*''2)} and 
SlHhim. All Cr3+ ions are equivalent in this spinel 
crystal structure and are expected to have the same mag- 
nitude, indicating that 5'i'^(i'0' 2 )) and 5'{''(2 > 2 ^0)} have 
to be collinear and orthogonal to each other. Among the 
spin structures shown in Fig. 6 and 7, the only possibil- 
ity would be the combination of S'{''(i''3'2)} shown in Fig. 
6 (c) and S'^'^^s-j-O)} shown in Fig. 7 (a). Fig. 8 shows 
the resulting coplanar and noncoUinear spin structure, in 
which each tetrahedra has two pairs of antiparallel spins 
and have zero net moment. 

For the k = (i, i, i) domain, there are many possible 
spin structures. From Figure 3 (b)-(d), we see that the 
intensities of the k = (i, i, i) reflections relative to other 
propagation vectors varies from sample to sample. These 
reflections are strongest for sample 1 shown in Figure 



FIG. 7: Prototypes of possible spin structures in a chemical 
unit cell for k = (i,i,0) in which satisfy the constraints 
described in the text: (a) -|-#''2ri ^ (1 + j)^kir2 ^ - 



k2T2 



klT4 _ 
/,kiT2_ 



{l-i)i>''^^^ -irP^^^i -^P^^^i , (c) iV"^^! +rp^^^^ +{l-i)i> 
{l + i)iP^^^^ - V''^^2 - iip^^^-^ + Tp^^^'^ + irp^^^'^ +{l + i)i>^^^^ + 
(1 - ■lytp^^^'' - iiP^^'^'i - ip^^'^'*, and (d) ii>^^^'i + ip^^'^'^ + {I - 

i)tp^^''^ + {1 + i)tp^^'^^ - lp^^'^'2 - iip^^'^'2 - itp^^-^'s _ ^karg ^ _ 

i)^''^"'' + (1 + i)ip^^^'' - V'*'"^* - Here ki = (i, i,0) 

and k2 = (| , ^ , 0) . The magnetic unit cell is doubled along 
the ar and the b-axes, and spins change sign in the chemical 
unit cell displaced by (1,0,0) or (0,1,0). Symbols representing 
z-coordinates of the magnetic Cr^"*" ions are the same as those 
in Fig. 6. 



3(b) and it is only for this sample that it is possible to 

distinguish between various models for the spin structure 
through goodness of fit. Thus, our discussion will be 
limited to sample 1. 

The basis vectors are given table VI. As the transi- 
tion is first order, multiple propagation vectors can con- 
tribute to the ordering. Many can be ruled out on phys- 
ical grounds, however, a large number of possibilities re- 
main. If one visualizes the B sublattice of the spinel 
lattice along the < 111 > direction, one can decompose 
it into alternating triangular and kagome layers. The ri 
configuration, has no spins in the triangular layer and 
the "q=0" spin configuration in the kagome layer. The 
T2 configuration has spins only in the triangular layer 
coupled antiferromagnetically. The T3 configuration is 
ferromagnetic in the kagome plane. The tq configuration 
allows for spins in the triangular layer to lie within that 
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FIG. 8: A combination sf^CLO.!)} + 5'{k(i,i,o)} j-^^. 
ZnCr2 04 which is consistent with all the experimental data 
explained in the text. Here S^^'-^'°-i^^ = V''^^" + tp'^^''^^ + 

i^kiTji _ -^kaTj^ _ j^kirsi _ j^k2T22 _ jfM^il + .(^'^2^22 with 

ki = (1,0, i) and k2 = (|, 0,1,0) (shown in Fig. 5 (c)), and 

5.{k{i,i,0)} ^^kir; _^ .^k^r; .^^kirs .^^k^ra _ 

i^kir^ _ ^kar^ ^.^jj ki = (i, i,0) and k2 = (1, 1,0) (shown 
in Fig. 6 (a)). 



plane at an arbitrary angle. We also considered linear 
combinations of these configurations. 

In table VII, we show the relative goodness of fits for 
the various spin models. Overall, we found that the 
model most consistent with our data was one with the 
"q=0" configuration in the kagome layer and an antifer- 
romagnetic configuration within the triangular layer (see 
Fig. 9). Though the difference in the values of is 
slight, if we examine a nuclear and a magnetic peak for 
goodness of fit for two different configurations, ri + T2 
and Ti + T2 + T3 models, we can see that there is is a real 
improvement in the fit for the ti + T2 spin model shown 
in Fig. 9 (see Fig. 10 and Table VII). 



2.0 



L5 



1.0 



0.5 



0.0 I 
0.0 







^ 1 I4J 1 1 ■ 




X 


^\ ^ 

j3 \ 

i^^ E 




M^ 


^ M 

M^ c 


^\ 


\ 




^ M 

c 


^ M 
M^ 


M 


^ \ 

1^^ E 


^ M 




^ M 


B \| 










C 



0.5 



1.0 



1.5 



3.0 



FIG. 9: This is the magnetic structure most consistent with 
the data for the k — (^, ^, 5) propagation vector. The ver- 
tices represent spins. Spins at vertices A,B, and C lie in the 
kagome layer and lie within the plane and have basis vectors 
of (0,1,1), (1,0,1), and (1,1,0) respectively. The vertices 
labeled by S lie in the triangular layer and the spins point 
along the body diagonal (1,1,1). The spins in the triangular 
layer form an antiferromagnetic pattern, alternating pattern 
out/into of the < 1, 1, 1 > plane. 



D. Summary 



We have solved the magnetic structure of ZnCr204. 
We have examined single crystals and three polycrys- 
talline samples. From this, we have found that the system 
has different domains formed from different k- vectors. 
The relative phase fractions vary from sample to sample 
(sec Table IX). The k = (1,0,0) is the same collincar 
structure as that of ZnV204- The k = (1,0,^) and 
bfk = (|, ^, 0) domain has the in plane structure found 
in Fig. 8. The structure of the k = (i, i, 5) domain is 
found in Fig. 9. Since the antiferromagnetic transition 
is first order, multiple characteristic wave- vectors are al- 
lowed unlike in the usual case of second order magnetic 
transitions. This multiplicity of domains with different 
characteristic wave vectors suggests that even the ordered 
state of this frustrated magnet is degenerate. 



IV. MAGNETIC CORRELATIONS IN 

ZnCr2-2xGa2x04 



In this section, we study how magnetic correlations 
change with nonmagnetic doping. 
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FIG. 11: Contour map of inelastic neutron scattering for Q = 
1.5 A^^. (b) T— dependence of magnetic Bragg scattering 
from a powder (squares), am ~ J I{Q,uj)ATvQ^dQdT'iLu 

where Vm is the volume per Cr^"^ ion, and of lattice strain 
along a and c (circles) measured by single crystal neutron 
diffraction. The figure is reproduced from Ref. []]]. 



B. Second order transition to Neel state in 

ZnCr1.9Gao.1O4 



FIG. 10: The peak to the left is the k = (|, |, |) magnetic 
reflection. The peak to the right is the nuclear (1, 1, 1) peak 
for comparison, (a) the flt to the spin configuration formed 
from a linear combination of ri and T2. (b) fit to the configu- 
ration formed formed from a linear combination of ti,T2, and 
T-3- 



A. First order transition to Neel state in ZnCr2 04 



For completeness, we start with the phase transition 
in pure ZnCr204|l|. Fig. 10 shows that in the pure 
ZnCr204 long range antiferromagnetic order (squares 
in frame (b)) and the local spin resonance (frame (a)) 
appear simultaneously in a spectacular first order phase 
transition. It also shows that the magnetic order is ac- 
companied by a cubic to tetragonal lattice distortion (cir- 
cles in frame (b)). The tetragonal distortion lifts some of 
degeneracy due to geometrical frustration and allows the 
system to order magnetically. Furthermore the ordered 
state pushes spectral weight in the energy spectrum up to 
the local spin resonance at fuv « 4.5 meV. It is unusual 
that a long range ordered phase can support a local spin 
resonance. 



The weak nonmagnetic doping in ZnCr1.9Gao.1O4 does 
not change the nature of the low T phase. As shown in 
Fig. 11, below Tn ~ 12.5 K magnetic long range order 
(squares in Fig. 11 (b)), tetragonal distortion (circles in 
Fig. 11 (b)) occur along with the appearance of the local 
spin resonance at Tiu! « 4.5 meV (Fig. 11 (a)). Mag- 
netic peaks in the doped material are the same as those 
in ZnCr204, which indicates that 5% doping of nonmag- 
netic Ga'^+ ions into Cr sites does not change the spin 
structure in the ordered phase. However, the three fea- 
tures appear gradually in a second order fashion, which 
is in contrast with the first order phase transition in the 
pure ZnCr204 shown in Fig. 10. Fig. 11 (c) shows 
that as soon as the static moment develops FWHM of 
magnetic peaks becomes Q'^esolution limited. This in- 
dicates the static correlations are long range no matter 
how small the static moment is. We conclude that the 
magnetic ordering in the weakly doped material immedi- 
ately develops in the entire material rather than in small 
magnetic clusters, growing in size grows as T decreases. 

Fig. 11 (a) shows that the local spin resonance and 
the low energy cooperative paramagnetic spin fluctua- 
tions coexist over the temperature range, 10 K < T < 
Tjv =12.5 K (also see Fig. 12 (b)). For comparison, in 
the pure ZnCr204 the Neel phase has only the linear 
spin waves below the local spin resonance (see Fig. 10 
(a)). To quantitatively study how the dynamic spin fluc- 
tuations in the weakly doped system evolve with T, we 
have fit the of / (Q=1.5 A.~^,hu}) in Fig. 11 (a) to two 
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FIG. 12: Contour maps of inelastic neutron scattering 
measured at Q = 1.5 A^^ as a function of energy trans- 
fer Twj and T. (b) T-dependence of the ordered moment 
(filled squares) and of crystal strains (opend and filled cir- 
cles) obtained by Rietveld analysis neutron powder diffraction 
data taken at BTl, NIST at various Ts using GSAS (c) T- 
dependence of fuU-width-of-Half-Maximum (FWHM) of the 
magnetic (1/2,1/2,2) Bragg peak (squares). The line is in- 
strumental angular resolution. 



simple non-resonant response functions, each with single 
relaxation rate: one centered at ?iw = with a relaxation 
rate Fi and the other centered at TujJq k, 4.5 meV with a 
relaxation rate F. 



i{nuj) = 



?i^0)2 + (|)2' 



(8) 



The first term is to account for the quasi-elastic scatter- 
ing that exists at T > 10 K. Since we did not have data 
above 30 K, it is difficult to extract meaningful infor- 
mation on the low energy excitations due to cooperative 
paramagnetism. Here we focuse on the local spin res- 
onance. Fig. 12 (d) - (f) show the results of the fits. 
The peak position, Tiui w 12.5 meV, (Fig. 12 (e)) is 
T-independent below T/v within experimental accuracy. 
The relaxation rate of the local resonance decreases as 
T decease to F = 1.4(1) meV at 1.4 K. For comparison, 
in pure ZnCr204 F = 1.5(1) meV for all T < Tn- The 
strength of the local resonance, Iq (Fig. 12 (d)), develops 
as proportional to the staggered magnetization, < M >, 
shown is Fig. 11 (b). This suggests that the static spin 
component is necessary to support the local spin reso- 
nance. 

Fig. 13 shows the spatial correlations of the fluctuating 
moments with different energies. The low energy lying 
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FIG. 13: (a)-(c) ?it<;-dependence of the inelastic magnetic 
neutron scattering intensity measured at Q — 1.5A^^ at three 
different Ts spanning the phase transition. Solid lines are the 
fits described in the text, (d) T-dependence of integrated 
intensity in the unit of (1/meV/Cr), (e) of peak position, 
and (f) of fuU-width-at-half-maximum of the fiu! m 4.5 meV 
excitations. 
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FIG. 14: Q-dependence of the inelastic magnetic neutron 
scattering intensity at 1.7 K and 12 K which is integrated 
over different energies. 



excitations at 12 K (squares) and the local resonance at 
1.7 K (filled circles) and 12 K (open circles) have almost 
identical Q-dependence with a Half Width at Half Max- 
imum K = 0.50(5) = 0.67(6)o*. Even though their 
characteristic energies arc different, the structure factor 
associated with the spin fluctuations have the same wave 
vector dependence. The excitations for < 3 meV 
at 1.7 K (triangles) also have a broad peak centered at 
Q =1.5 A-^ 
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C. Spin freezing in ZnCri,4Gao.604 



neutron scattering data 



0.3 



I" 




l9' 



/a2/-i iHQ)/\FiQ)WdQ 



<M> ' : . (10) 



Integrating the difference data over Q yields | < Af > 
|2 = 0.24(4)^|/Cr, in other words, | < M > | = 5 < 
S > fiB = 0.49(4)(UB/Cr. This quantity is substantiaUy 
less than the Neel value | < M > | = gS^B — S/xs/Cr 
and also much less than those values of ZnCr204 and 
ZnCr1.9Gao.1O4. 
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FIG. 15; Q-dependence of elastic magnetic neutron scattering 
intensity measured with energy window of \hui\ < 0.05 meV. 



In this section, we study spin correlations in 
ZnCr1.4Gao.6O4 which exhibits spin-glass-like behaviors 
in bulk susceptibility measurements |4(l|. Fig. 14 shows 
elastic magnetic scattering intensities measured at 1.4 K. 
High temperature background was measured at 20 K and 
subtracted. Unlike in the pure ZnCr204 and the weakly 
doped ZnCr1.9Gao.1O4, this system does not have mag- 
netic Bragg peaks but a broad peak centered at a finite 
wave vector Q = 1.5 A"^ with k = 0.48(5) This 
indicates that the 30% nonmagnetic doping destroys the 
magnetic long range order and reduces the correlation 
length down to distance between nearest neighboring Cr 
ions. I{Q) going to zero as Q approaches zero indicates 
that the antiferromagnetic constraints are still satisfied 
in the heavily doped sample. The solid line is the fit to 
the powder-averaged magnetic neutron scattering inten- 
sity for an isolated spin dimer[53ll55|. 



HQ) « \F{Q)\ — q;^^ — 



(9) 



where the distance between nearest neighboring Cr"^+ 
10ns ro = 2.939 A. The spin pair model produces a 
broader peak than the experimental data. Instead, the 
Q-dependence is almost identical to that of the fluc- 
tuating spins in pure ZnCr204 and weakly diluted 
ZnCr1.9Gao.1O4 (see Fig. 13). This suggests that the 
same local spin objects involving more than isolated 
dimers are responsible for the broad Q-dependence in 
those materials even though the energy for the correla- 
tions changes with doping of the nonmagnetic ions. We 
can estimate the average frozen moment from the elastic 
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FIG. 16: T-dependence of elastic magnetic neutron scat- 
tering intensity interaged over < 0.05 meV and inelastic 
magnetic neutron scattering intensity integrated over 0.3 meV 
< huj < 0.6 meV. Both data were measured at Q = 1.5 A~^. 



Fig. 15 shows elastic and inelastic neutron scattering 
intensities measured at Q = 1.5 A~^ as a function of 
temperature. The onset of elastic magnetic scattering at 
around 10 K signals the development of magnetic cor- 
relations on a time scale, t > 2h//S.E — 0.013 ns set 
by the energy resolution of the instrument. Energy inte- 
grated inelastic scattering data over 0.3 meV < huj < 0.6 
meV show a broad maximum at around 10 K indicating 
the critical scattering at the phase transition. Bulk sus- 
ceptibility data with a maximum at a lower temperature 
Tf — 1.8 k|4^ show that this apparent critical tempera- 
ture is not unique but decreases with the energy scale of 
the measurement. Such behavior, which is characteris- 
tic of spin glasses, indicates that a precipitous softening 
of the magnetic fluctuation spectrum takes place upon 
lowing the temperature, leading to anomalies when the 
lowest energy scale of the system falls below the charac- 
teristic energy scale of the measurement. 
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V. DISCUSSIONS AND SUMMARY 

ZnCr204 is so far the best realization of an antiferro- 
magnet on the magnetic lattice of corner-sharing tetra- 
hedra with uniform nearest neighbor interactions and 
without any site-disorder. Upon coohng, this system is 
heading toward spin-hquid state with the signature of 
almost linear spin relaxation rate. At low temperature 
T < Ttv = 12.5 K , the system undergoes a cubic-to- 
tetragonal distortion to settle into a Neel phase with a 
local spin resonance [l|. The three features, tetragonal 
distortion, long range order, and the local spin resonance, 
occurs abruptly in a first order fashion. The spins in 
the Neel phase have reduced staggered magnetization, 
< 1 (see the inset of Fig. 2) due to geometrical 
frustration. Weak 5% nonmagnetic doping into the mag- 
netic lattice further suppresses the staggered magnetiza- 
tion but does not destroy the Neel phase at low tem- 
peratures. The phase transition from cooperative para- 
magnetic phase to Neel phase, however, occurs gradu- 
ally upon cooling in a second order fashion. The cubic 
to tetragonal lattice distortion also follows the develop- 
ment of the magnetic phase transition. This supports 
that the phase transition is magnetically driven. It is 
understandable that the nature of the low temperature 
phase does not change with 5% doping because for 5% 
dilution in the magnetic lattice, the majority of tetrahe- 
dra have all 4 spins (81% of tetrahedra have all 4 spins 
and 17% have 3 spins[P|). For 30% dilution (l-x=0.7), 
24% of tetrahedra have 4 spins, 41% have 3 spins and 
27% have 2 spins. Even though it is still above the per- 
colation threshold {1 — Xc — 0.39), the long range corre- 
lations are destroyed and replaced with short range ocr- 
relations. Despite these differences, all three materials 
contain spin correlations with a common broad Q depen- 
dence even though the energetics of the local correlations 
change with the occupance of the magnetic lattice and 
the existence of long range order. This indicates that 



the local spin object responsible for the common Q de- 
pendence is robust against strong disorder. This finding 
may explain why bulk properties in geometrically frus- 
trated magnets are robust against dilution. In SCGO(x), 
where the magnetic entity relevant to geometrical frus- 
tration can be viewed as quasi-two dimensional (111)- 
slabs of corner-sharing tetrahedrafssj, the bulk suscepti- 
bility shows field hysteresis and the nonlinear susceptibil- 
ity diverges, typical of spin glasses, the specific heat C{T) 
is proportional to as in an ordinary two-dimensional 
antifcrromagnct[56j. These bulk behaviors are very ro- 
bust against magnetic dilution. [56], |53| Our finding in- 
dicates that the fundamental spin degree of freedom in 
ZnCr2-2xGa2x04 is the hexagonal loop of antif erromaq- 
netic spins observed in the pure ZnCr204.'24'| The low 
energy physics is governed by the excitations of the local 
spin degree of freedom and therefore is robust to dilution. 

In summary, we have determined the spin structure of 
the Neel phase in ZnCr2_2xGa2x04 which would pro- 
vide a starting point for a theory for this system. We 
have studied, using neutron scattering, how nonmagnetic 
doping changes the first order magnetoelastic phase tran- 
sition in pure ZnCr204 into the second order spin-glass- 
like phase transition. We have found that a broad Q 
dependence is robust against dilution, suggesting that 
such local spin correlations both in Neel phase and in 
short range ordered phase is intrinsic to the geometri- 
cally frustrated magnets and distinguishes these systems 
from the ordinary spin glasses. 
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TABLE I: Measured ratios of SF to NSF scattering intensities for three magnetic reflections that belong to k = (1,0, |' 

{h,l,k) (crSF/o'NSF)obs ^ — 

(1,0.5,0) 0.4(1) 0.2 
(1,1.5,0) 0.7(1) 0.6923 
(1,2.5,0) 0.83(5) 0.862 
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TABLE II: The basis functions of the irreducible representation of group Oj^ for the two-arm star of the wave vector k = (5,^,0). 
Here k is in terms of the cubic notation and k = j{2, 1, 1) in Kovalev's notation [s^]. Here ki = k and k2 — — k. The notation 

of representations, such as ti, ti and so on, followed the Kovalev's notation, ipf^i'^ = li 2,3,4) represent the basis functions for 
the spins located at (0.5,0.5,0.5), (0.5,0.75,0.75), (0.75,0.5,0.75), and (0.75,0.75,0.5) in the cubic notation, respectively. This 
table was taken from Table 22 in page 131 of Ref. 50]. 
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TABLE III: The basis functions of the irreducible representation of group Fd3m{Ol) for the two -arm star of the wave vector 
k = (1,0, I). Here k is in terms of the cubic notation and k = j(l, 1,0) + |(0, 1, 1) in Kovalev's notation using primitive 
reciprocal unit vectors [H^l- Here ki = k and k2 = — k. The notation of representations, such as ri, Tj and so on, foUowed the 
Kovalev's notation, ipf'^i'i = 1)2,3,4) represent the basis functions for four spins in the primitive ceU located at (0.5,0.5,0.5), 
(0.5,0.75,0.75), (0.75,0.5,0.75), and (0.75,0.75,0.5) in the cubic notation, respectively. 
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TABLE IV: Superpositions of two irreducible representations, Citp + C2tp which give real spins to the atoms for ki = 
(5, |, 0) and k2 = — ki. 1(1, 1, 0) represents the spin at site 1 is along (1, 1, 0). Positions of the 1 to 16 sites are shown in Fig. 
6 (a). 

Superposition Nonzero spins in a chemical unit cell 
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TABLE V: Superpositions of two irreducible representations, Citp ^'^ + Citp which give real spins to the atoms for ki 
(1,0, i) and ka = -ki. 
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TABLE VI: The basis functions of the irreducible represen- 
tation of group Fd3m{Ol} for the two-arm star of the wave 
vector k = (i, i, i). 
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T5 


(0, 0, 0) 

(0,0,0) 
(0,0,0) 


(2,0,0) 

(0,2,0) 
(0,0,2) 


(0,1,0) 

(0,0,1) 
(1,0,0) 


(0,0,1 

(1,0,0 
(0,1,0 




(0,0,0) 

(0,0,0) 
(0,0,0) 


(1,0,0) 

(0,0,1) 
(0,1,0) 


(0,2,0) 

(2,0,0) 
(0,0,2) 


(0,0,1 

(0,1,0 
(1,0,0 




(0,0,0) 
(0, 0, 0) 
(0,0,0) 


(1,0,0) 
(0,1,0) 

(0,0,1) 


(0,1,0) 
(0,0, 1) 

(1,0,0) 


(0,0,2 
(2,0,0 

(0,2,0 


T6 


(2,1,1) 
(0,1,1) 


(0,0,0) 
(0,0,0) 


(0,0,0) 
(0, 0, 0) 


(0,0,0 
(0, 0, 
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TABLE VII: Goodness of the fit of the different spin models 
with k = (|, |, ^) to the , ^, ^} reflections obtained from 
different ZnCr2 04 poly crystalline samples. The best fit was 
obtained with the ti + T2 model (see Fig. 10 (a)). 



Sample 


Rep. 


T (K) 


x' 


Rwp 


1 


Tl 


1.5 


4.915 


.0664 




Tl,T2 




4.888 


.0662 




Tl,T2,T3 




4.911 


.0664 




Tl,T6 




4.998 


.067 




T2,T3 




4.995 


.067 


2 


Tl,T2 


7.3 


1.368 


.0618 


3 


Tl,T2 


2 


4.003 


.0592 



